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Abstract 

Recently, Mueller and Son discussed the time evolution of a dense 
system towards equilibrium in a scalar Xip^ field theory PJ . They show 
the equivalence of the classical field approximation and the Boltzmann 
equation in all but linear terms in the occupation number. Here we 
present the generalization to high energy QCD. 
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1 Introduction 



In high energy heavy ion colhsion experiments it is beheved that a new state 
of matter, the so-called quark gluon plasma, may be created (j2], Theo- 
retically it is therefore important to understand the time evolution of dense 
strongly interacting matter towards thermal equilibrium. This evolution can 
be described in terms of the gluon occupation number f. In the very early 
stages after a collision f is of the order a^"^ (j3j, jS]) where as is the strong 
coupling constant. As the system evolves f decreases and is finally of the 
order 1 at equilibrium. 

Unfortunately there is no general theoretical framework to describe the 
complete evolution. For large f (i.e. very early times) non-linear classical 
field theory can be applied while a Boltzmann equation that takes quantum 
fluctuations into account (this will be called full or "quantum" Boltzmann 
equation) may be used near equilibrium. 

In order to obtain a continuous description of the evolution one may try 
to fit the transition between the two approaches by a classical field approx- 
imation of the full Boltzmann equation (Fig. 1). However, it is a priori not 
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Figure 1: time evolution of the gluon occupation number 

clear that these transitions are as smooth as the above figure suggests. In 
this paper we compare the classical and quantum Boltzmann equations in 
high energy QCD and show that they are equivalent to orders and 
The comparison of non-linear classical fields with the classical Boltzmann 
equation is an entirely different problem and will not be treated here. 
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The structure of the paper is as follows: in section 2 we briefly review the 
main results of Mueller and Son where the basis for our considerations 
is developed within a scalar field theory. In section 3 we generalize their 
strategy to high energy QCD and show that one obtains essentially the same 
results. 

We want to emphasize that our argumentations are mostly qualitative. 
Perhaps surprisingly, we need very few computations to derive our results; 
some simple topological and diagrammatical considerations are sufficient. 
Similar methods can be found, e.g., in p. 



2 Scalar field theory 

2.1 Separation of classical and quantum contributions 

Our starting point is the scalar field theory with X^p'^ interaction given by 
the Lagrangian 

C = \d,^d^^-\m'^'-^^' . (1) 

Throughout this paper we assume coupling constants to be small and perform 
the relevant computations in first order perturbation theory. 

Our ultimate interest will be to gain deeper insight into the time evo- 
lution of heavy ion collisions. So we have to generalize our theory to finite 
temperature. It is convenient to apply the closed time path formalism (OTP) 
which leads to a doubling of the field variables (jZj, (Hj): 

The Lagrangian then reads 

c-cTP = ^a^$-<9^$- - ^m2$2 - A$4 

2 2 4! j^2) 

We now perform a change of the field variables in order to distinguish the 
classical field (denoted $) and quantum fiuctuations (denoted 11) 

$ = i($- + $+) n = <l>_-<l'+ , (3) 

$_ = $ + In $+ = $ - . (4) 

As we are interested in systems where both $+ and $_ are large, the above 
interpretations of $ and 11 are already at least quahtatively justified. 
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In terms of these new fields the Lagrangian becomes 

= a^^a'^n - m^^n - ^('^'^n + ^^n^) . (5) 

Indeed, it can be easily shown now by neglecting higher than linear terms in 
n that $ fulfills the classical equation of motion (for details see P) 

(□ + m2)<l> + ^<l> = . (6) 
From the Lagrangian (5) one can identify the vertices of the theory (Fig. 2). 



Figure 2: vertices of the scalar theory 



The full lines correspond to the classical field $ and the dashed lines to the 
quantum fluctuations 11. The diagram on the left will be called classical 
vertex, the right one quantum vertex. 

The most important quantities for our argumentations will be the free 
Greens functions or propagators for the $ and 11 fields, respectively. In a 
more rigorous discussion one should adopt the full propagators, but mass 
corrections are negligible as long as the coupling constant is small enough, 
i.e. /A <^ 1. (For details on this point see [T].) 

Note that there are also mixed propagators describing the change from $ 
to n and vice versa. One obtains 

G^^ = 2n6ip'-m')if + ^) , (7) 



G 



u<s> 



p2 _ ^2 _ j^^p^ 



G<s>u — -5 5—- — ) (9) 

G'nn = . (10) 

The crucial observation here is that only 6*$$ depends on the occupation 
number f. 
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2.2 Boltzmann equation 

The Boltzmann equation describes the time evolution of the occupation num- 
ber f of a given state as the difference between the scattering of particles into 
and out of this state (gain and loss). These scatterings are included in the 
collision term C. Considering the scalar field theory (1) one may express the 
collision term to lowest order diagrammatically (Fig. 3). 




"gain" "loss" 



Figure 3: collision term to lowest order 

For a detailed analysis of the Boltzmann equation in scalar and gauge field 
theories see j^]. For our purposes it is sufficient to examine the topology of 
the involved diagrams and it is easy to see that these are effectively sunset 
graphs (Fig. 4). 




Figure 4: collision term topology 



When we express now our theory in terms of $ and 11 we use the corre- 
sponding vertices. The classical approximation then consists in retaining 
only collision term diagrams without quantum vertices. 
There are three different classical diagrams (Fig. 5). 




Figure 5: classical contribution to C{p) 



In Fig. 5 the left diagram is proportional to (/+|)^ as three propagators 
appear. (Here and in the following empty circles appear in propagators while 
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full circles stand for vertices.) Similarly, the other graphs are proportional 
to (/ + and (/ + ^Y, respectively. 

In the full Boltzmann equation one also has to include the diagrams with 
one quantum vertex (Fig. 6). In Fig.6 the left diagram is proportional to 




Figure 6: quantum contribution to C{p) 



(/ + while the others are independent of f. 

There are no first order diagrams with two quantum vertices as at least 
one propagator Gnn = would appear. So we can conclude that classical 
field theory and the Boltzmann equation are equivalent in orders and 
(PJ). In the next section we will show how this strategy can be generalized 
to QCD of heavy ion collisions. 

3 High energy QCD 

3.1 Separation of classical and quantum contributions 

The results of the previous section may be generalized to QCD by making 
two approximations which are well established in the context of heavy ion 
collisions. 

First, we will work within the gluon saturation scenario ([IHl, [12] i 
|13j). This means that we may neglect the fermionic degrees of freedom in 
our system. So the QCD Lagrangian simplifies considerably and reads in the 
CTP formalism 

CcTP = -\FrF,,Ma,] + \FrFa,AK'\ ■ (11) 

Here the letter 'a' is a color index and the gluon field. In analogy to the 
field transformation in the scalar case we now define 

^a/j = -^{^a^i + = A^^ — A^^ , (12) 

^a/. = -^-aM + ^n,^ /l+=<l>,^-^n,^ . (13) 



5 



This leads to the following expression for 



(14) 



and similarly for 

Next we express the Lagrangian (11) in terms of the fields $^ and 11^. 
Neglecting higher than linear terms in the quantum fluctuations one ob- 
tains 



Ainear ^ (i^^^^ri^Dn^ (15) 

with 

J^a^^Q^^a, _^^a,c^c^ , (16) 

F^^'l^ = di'^ - + gf"'^^"^^"^ . (17) 
Thus, fulfills the classical equation of motion 

<^r[^]=o . (18) 

So our interpretation of as classical field is justified in complete analogy 
to the scalar case. 

Prom the equations (15) to (17) one also sees that there is only one 
classically allowed vertex in first order in the coupling constant g (Fig. 7). 




Figure 7: classical vertex in QCD 



Here full lines denote classical fields and dashed lines quantum corrections 
like in the previous section. 
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In order to obtain the quantum couplings one has to take into account 
the terms nonhnear in 11^. The corresponding Lagrangian may be easily 
computed and reads 

>C„onz.„ = -IgfatcTl'^U^" ilid.K - dull';) + gr'^^iK) ■ (19) 
So we have one first order quantum coupling (Fig. 8). 

< 



Figure 8: quantum vertex in QCD 



3.2 The collision term 

Our second approximation leads to a simple topology of the collision term 
C. As fermions are neglected the only contributions to C come from gluon 
scattering. We now assume that the t-channel dominates, where t is the 
Mandelstam variable (see, e.g., |^). (In fact, our only assumption is to 
work in the high energy limit where both approximations are valid.) 
The t-channel gluon scattering is shown diagrammatically in Fig. 9. 




Figure 9: t-channel gluon scattering 

This makes the relevant collision term topology quite simple as can be seen in 
Fig. 10 (note that in contrast to the previous section only three-field vertices 
are present). 

As the main features of the propagators do not change, especially 6*$$ oc 
(/+|) ^^'^ Gun = 0, we are now again ready to compare classical field theory 
with the full Boltzmann equation. As before the classical approximation 
consists in neglecting diagrams with quantum vertices which are only taken 
into account in the full collision term. 
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Figure 10: collision term topology in QCD 



Let us consider for example the diagram in Fig. 11 that clearly has the 
required topology 




Figure 11: quantum contribution to C{p) 

It contains one quantum vertex so it is not included in the classical approxi- 
mation. The contribution from this graph to the collision term is proportional 
to (/+|) as one propagator 6*$$ appears. Similarly, the diagram in Fig. 12 is 
classical and proportional to (/ + ^)^- It is easy to check that these examples 




Figure 12: classical contribution to C{p) 

give already the highest order contributions to the classical approximation 
and the quantum corrections, respectively. 

We conclude that the results of the scalar case are indeed valid in high 
energy QCD: The classical field approximation and the quantum Boltzmann 
equation match in orders and i.e. in leading orders of the gluon 
occupation number. 
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